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A two-dimensional gas of non-interacting quasiparticles in 
a nearly periodic potential is considered at zero temperature. 
The potential is a superposition of a periodic potential, in- 
duced by the charge density wave of a Wigner crystal, and 
a weak random potential due to disorder. There is a metal- 
insulator transition that is controlled by the strength of the 
periodic potential. The transition is continuous in the pres- 
ence of randomness. We evaluate the density of states, which 
is non-zero at the Fermi energy in the metallic phase, and the 
dc conductivity. The latter changes with decreasing modula- 
tion of the periodic potential from to a « 2e 2 /h. 

PACS numbers: 71.30. +h, 73.40.Qv 

The metal-insulator transition in silicon MOSFETs Jl]] 
and GaAs/AlAs heterostructures presents a new, 
unexpected and interesting phenomenon in a two- 
dimensional electron (or hole) gas |l| . One of the charac- 
teristic features of the transition is the scaling behavior 
of the resistivity with temperature. There is a critical re- 
sistivity which shows up at the transition point from the 
insulating to the conducting regime. The corresponding 
critical conductivity typically ranges from 0.3e 2 //i |l| to 
2e 2 /h j|. In other words, it seems from the experiments 
that at zero temperature there is either an insulating 
state or a conducting state with a conductivity equal or 
larger than a critical (minimal) conductivity. The na- 
ture of the conducting state (high electron density) is 
unclear. It is either a normal metal a superconduct- 
ing state [Q or a state controlled by charged traps [^). 
The nature of the insulating state (low electron density), 
on the other hand, is less controversal. Since the 2D elec- 
tron gas must have a high mobility in order to undergo a 
transition to a non-insulating state, disorder is very weak 
in the samples. Therefore, a Wigner crystal or at least 
a modulated electron density with short range order is 
expected due to the Coulomb interaction ||. In the fol- 
lowing we shall discuss a model which provides a simple 
picture of a metal-insulator transition in a nearly peri- 
odic potential created by a modulated electron density. 
The discussion will focus on the approach of the conduct- 
ing regime, coming from the insulating regime, but does 
not include a description of this regime away from the 
transition point. 

Lattice dynamics calculations [0J|] and computer sim- 
ulations for the pure 2D electron gas indicate that 
the Wigner crystal forms a hexagonal lattice. Quasipar- 
ticles, which are the excitations in the Wigner crystal, 
experience an effective potential due to the modulation of 



the electron density. Formally, the dynamics of the quasi- 
particles can be derived from the model of a 2D electron 
gas which is subject to Coulomb interaction. Starting 
with a microscopic model we can apply a self-consistent 
approximation for the space-dependent electron density 
n (cf. |ll]]), which describes the Wigner crystal, and re- 
gard the fluctuations around the static electron density 
as quasiparticles. In leading order the interaction of the 
fluctuations is neglected, i.e. we consider independent 
quasiparticles. This approach is analogous to the deriva- 
tion of the Bogoliubov de Gennes Hamiltonian in a su- 
perconductor, where the superconducting order parame- 
ter A is treated in a self-consistent (BCS) approximation. 
In contrast to A the electron density n couples directly 
to the quasiparticle density. Therefore, it can be consid- 
ered as an effective potential for the quasiparticles. The 
latter prefer to stay in places with low electron density 
(minimal potential, the circles shown in Fig. 1.) but can 
also tunnel through the saddle points of the potential 
between these potential minima (dashed lines in Fig. 1). 




FIG. 1. Effective potential and the lattice of the 
tight-binding model. Circles indicate minima of the electron 
density and dashed lines the tunneling of the quasiparticles. 

The model of the quasiparticles can be further sim- 
plified by assuming a discrete potential instead of the 
smoothly varying density n. In other words, the vari- 
ation of n is taken into account only down to a finite 
length scale. For the latter we choose the distance be- 
tween a minimum and an adjacent saddle point. This is 
a minimal model in order to study the broken transla- 
tional invariance caused by the density modulation. It 
can be considered as a tight-binding approximation for 
the quasiparticles, where the lattice points are the min- 
ima and the saddle points of the electron density. There 
is a potential Vr at each lattice site R with Vr = V+ 
on the minima and Vr, = V- on the saddle points. For 
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the sake of simplicity of our model, (I) only the overlap 
between nearest neighboring lattice sites (i.e. a mini- 
mum and a saddle point) and (II) no spin effect are 
taken into account. It is convenient to substract a con- 
stant potential in order to replace Vr by the modulation 
tor = Vr — (V+ + V_)/2. Then tor is a staggered field 
on the lattice with tor = ±(V+ — V-)/2 = ±m. This 
field breaks the translational invariance on the lattice. 
A translational-invariant representation can be obtained 
from one of the two sublattices. Using the sublattice of 
the maxima (circles in Fig. 1), the co-ordinates of a site 
on this sublattice are r = (x,y), where the x (y) direc- 
tion is horizontal (vertical) in Fig. 1. The co-ordinates of 
a site of the other sublattice (corners of the triangles in 
Fig. 1) are obtained by shifting the co-ordinates r with 
the unit vector a = (0, 1). Any site R on the lattice can 
now be written as (r,j) = r + (j — l)a, where j = 1, 2 is 
the index of the sublattice. 

The quasiparticle Hamiltonian is defined as 



H 



EE 

r,r' j,j'=l 



(1) 



where c (c') are the annihilation (creation) operators of 
the (fermionic) quasiparticles. For the off-diagonal ma- 
trix elements we assume H, 



r.r+e 



t, where the unit vec- 
tors e = a, b, c are indicated in Fig. 1. This Hamiltonian 
has also been studied in connection with the quantum 
Hall effect. 



The matrix H — Hq + ma^, where the 
Pauli matrix 173 refers to the sublattice index, can be 
given in Fourier representation with respect to the sub- 
lattices (r — > (k x , ky j) as 

H — > H = t[c(k x ,k y )ai + s{k x ,k y )a 2 ] + 77103 (2) 

with Pauli matrices a±, <j% and 

c(k x ,k y ) — cos(-ky) + cos(v / 3fc a; /2 + k y /2) 

+ cos (-V3k x /2 + k y /2). (3) 

s{k x , k y ) is obtained from this expression by replacing the 
cosine by sine. For this translational-invariant Hamilto- 
nian we obtain the dispersion 



E±(k x , k y ) = ±yjm 2 + t 2 {c 2 + s 2 ) 



(4) 



as shown in Fig. 2. The modulation to creates a gap 
between the two bands with E + > and with E- < 
0. This gap vanishes only in the limit of a vanishing 
modulation to = 0. Thus our system is insulating as 
long as we have a modulated electron density. In the 
case m = there are six nodes (k x , k y ) on a circle with 
radius 4 • 3~ 3 / 2 7r, located at 
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For the analysis of the transport (i.e. low energy) proper- 
ties of this system it is sufficient to study the properties 



associated with the nodes. Expansion around the two 
nodes at (±47r/3 3 / 2 , 0) leads to the Dirac Hamiltonians 



iJ ~±(3t/2)(A:i<7i±fc 2 o- 2 ), 



(6) 



and around the other four nodes to 

H - -(3*/4)[(V3fci + k 2 )ai + (V3ki - k 2 )a 2 ], (7) 

where k\ (k 2 ) is now the deviation from the correspond- 
ing node in x- (y-) direction. Thus for the pure system 
the Hamiltonian separates into the six nodes when we 
consider small energy (i.e. large scale) behavior, where 
the vicinity of each node is given by a two-dimensional 
Dirac operator. 




FIG. 2. Dispersion E(k x , k y ) of quasiparticles, measured in 
units of the overlap t. 

So far we have assumed that there is perfect long-range 
order in the Wigner crystal, leading to a translational- 
invariant Hamiltonian for the quasiparticles. This is an 
idealization which is not valid in a real system where we 
have always impurities. It is known that the 2D Wigner 
crystal can be strongly affected by impurities. In the 
presence of a high density of impurities the crystalline 
order is destroyed completely and an electron glass can 
be created On the other hand, short-range order 

may survive if disorder is only weak. This situation can 
be described in terms of the quasiparticles by a weak 
random perturbation added to the periodic modulation. 
Various types of disorder are possible but we will restrict 
the discussion to the simplest case, 



M r (T3 = (to + SM T )a 3 , 



(8) 



assuming that the effect of this type of randomness de- 
scribes a generic situation. 

In principle, the disorder SM can couple different 
nodes because the random term has fluctuations on all 
length scales. However, we will assume that the overlap 
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between the scales, related to the various nodes, is neg- 
ligible, i.e. for k and k! belonging to the vicinities of two 
different nodes we have 

J2e- ikT+ik ' r M T nO. (9) 

r 

For each node we introduce an independent random term 
5M a ^ T which means that M a>r fluctuates only on scales 
which are compatibel with the corresponding node. We 
assume a Gaussian distribution with (8M a ^ T ) — and 
(SM a<r SM a ty) = g8 a ,a'8 r .r'- Then the Green's function 
G(z) = (H — z)" 1 can be block-diagonalized with respect 
to the six nodes as diag[(ffj — (H e — z)^ 1 ) with 

independent Dirac Hamiltonians 

H a = z'(o a Vi + b a V 2 )cri + !'(c a Vi + g? q V 2 )ct2 + M a ^cr 3l 

(10) 

where a a = —d a = — 3i/2, b a — c a = for the nodes 
with k y = (i.e. a = 1,2) and a a = c a = — 3 3 / 2 </4, 
b a = —d a = — 3i/4 for the nodes with k y ^ (i.e. 
a = 3,..., 6). The block-diagonal structure is a cru- 
cial simplification because we only have to evaluate the 
physical quantities for the six nodes independently The 
density of states (DOS) and the conductivity of the cor- 
responding Dirac Hamiltonians are known and shall be 
discussed subsequently. 

The DOS of a system near a metal-insulator transition 
exhibits a characteristic behavior. In the simplest case 
there is a gap in the insulating phase which closes at the 
transition to the metal. The closing of the gap was mea- 
sured recently in a tunneling experiment near the metal- 
insulator transition of boron-doped silicon crystals, a 
three-dimensional electron gas. |14| ft was found that the 
DOS goes roughly like a power law p(E) ~ p${\E\ / E y 
(E is the distance from the Fermi energy Ep and 7 ss 0.5) 
in the insulating phase and like p(E) ~ p [l + {\E\lEo) 1 ] 
in the metallic phase. Thus the insulating phase can be 
distinguished from the metallic phase only for energies 
close to the Fermi energy. A similar behavior is expected 
for a two-dimensional system, with a different exponent 7 
though. For instance, a linear behavior p(E) ~ pq\E\/Eq 
was observed in a 2D electron gas subject to a perpendic- 
ular magnetic field. ]l5j This result is in agreement with 
the effective Dirac fermion description of quasiparticles 
in a quantum Hall system JlJ] which also gives a linear 
DOS in the absence of disorder. Moreover, a linear be- 
havior of the DOS was found in a mean-field calculation 
of a disordered 2D electron gas. |l3[ ] 

In our model, the block-diagonal structure enables us 
to write the DOS as a sum of DOS's of the vicinity of 
each node. Using the result of the Dirac fermions from a 
saddle point approximation this gives for m = at each 
node a [[l6) 

Pot{E )K e —^- + \E\\ (11) 



where the exponent 7 varies with increasing randomness 
g from 1 to (Fig. 3). 
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FIG. 3. Density of states at m = for different strength of 
randomness 

The energy scale is set by the overlap t. The DOS as 
a fu nction of the m odulation parameter m at E = is 
cx V4e~ 27r /s - m 2 ©(4e- 27r /f - m 2 )/2g. Therefore, the 
DOS vanishes with g — > 0. The effect of the periodic 
modulation is the creation of a gap. On the other hand, 
the random fluctuations around the periodic modulation 
create states at the Fermi energy as long as the random 
fluctuations are strong relative to the periodic modula- 
tion. To have quasiparticle states at the Fermi energy 
it requires random fluctuations, characterized by the pa- 
rameter 2e _1T / 9 , which are larger than the periodic mod- 
ulation, characterized by m, . The creation of new states 
around m = by the randomness can be explained by 
the formation of tail states at low energy. This is similar 
to the creation of Lifshitz tails. [M However, it will be 
discussed in the following that these states are conduct- 
ing - in contrast to the localized states in the Lifshitz 
tails. 

In the presence of a gap the conductivity vanishes with 
vanishing temperature. Here we are only interested in 
zero temperature properties, i.e. we must study the 
transport properties for the system with a closed gap. 
Since a randomly disturbed modulation (5M r ^= 0) cre- 
ates states at the Fermi level for m < 2e -7r / 9 , the case of 
the nearly periodic modulation is a candidate for a metal- 
insulator transition. However, the vanishing gap is nec- 
essary but not sufficient for a conducting state because 
the quasiparticle states might be localized by random- 
ness. The localization effect is known to be very efficient 
in d = 2, e.g., all quantum states are localized according 
to conventional scaling theory for localization [0 . How- 



3 



ever, it was found that Dirac fermions with random mass 
can escape from localization by a special mechanism due 
to symmetry breaking p8| . It leads to diffusion of the 
Dirac particle for sufficiently small average mass, i.e. for 
sufficiently small modulation. The conductivity can be 
evaluated using the Kubo-Greenwood formula. Then the 
dc conductivity at T = and E = reads 

a ps e - lim if V r 2 Tr 2 (G(0, r; i V )G(r, 0; -irf)) 

tl r)—>0 *• — ' 
r 

= -Li im y V £(7MU=o (12) 

ft rj— >0 

with the two-particle Green's function C(r,rj) = 
Tr 2 (G(0, r; irf)G(r, 0; — irf)). The latter is a sum of con- 
tributions from the six different nodes due to the block- 
diagonal form of the single-particle Green's function. The 
average two-particle Green's function can be diagonalized 
by a Fourier transformation because it is translational in- 
variant as a consequence of the uniform distribution of 
the randomness. The nodes require a similarity transfor- 
mation 

fci — > k[ = a a k\ + b a k2, fc 2 — * k' 2 = c a k\ + d Q fc 2 (13) 

in order to obtain the usual Dirac form a-k for the Hamil- 
tonian. The average two-particle Green's function is flg| ] 



C a (k,rj) 



D{k' 2 + k,' 2 2 )/\a a d a 



^a^a | 



(14) 



where D is the diffusion coefficient of the Dirac fermions. 
The extra factor l/\a a b a — c a d a \ is the Jacobian of the 
transformation k — > k' . For weak disorder (j « 1) the 
diffusion coefficient is 19 



D 



1 



2tt 2 p 



(1 



^)[1-^(1 

Z7T 



2x 2 )]e(l-x 2 ) (15) 



with x — me* 1 9 /2. With ( |l2"| ) and the summation over 
all nodes this implies for the conductivity 



cr(m) ps -r-s (l 
ft 



)[l-JL(l-2x- 2 )]Q{l 



(16) 



with s ~ 2.1. The conductivity depends on the modu- 
lation m which scales with e^l 9 . The transition is very 
sharp for weak disorder, since the scale is exponential, 
and it is discontinuous with a minimal conductivity 2e 2 /ft 
if disorder is absent. 

The quasiparticle states overlap in such a manner that 
they form conducting states. This effect is known for 
Dirac fermions with randomness already in one dimen- 
sion, where the zero energy modes are extended. ]2(i| ] 
However, in contrast to the finite diffusion coefficient 
( |l5| ) , leading to a diffusive behavior of quasiparticles, this 
coefficient is singular then. In our model the diffusion is 
valid only for weak randomness, as it is indicated by the 
vanishing diffusion coefficient (p5h if g = 2ir. 



In conclusion, a possible metallic state in a two- 
dimensional electron gas is approached from an insulat- 
ing regime by the destruction of the gap of the insu- 
lating state by disorder. Our discussion is based on a 
simple model for quasiparticles in a nearly periodically 
modulated electron density, which can be considered as 
a Wigner solid with weak disorder. The system becomes 
metallic for a sufficiently weak modulation with a con- 
ductivity a ps 2e 2 /ft. The latter decreases linearly with 
increasing disorder. The conductivity varies continuously 
for the nearly periodic potential. 
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